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ABSTRACT 

In previous studies, it has been shown that the large line widths observed in high surface brightness Ho- emitters at low and high 
redshifts are likely due to the mechanical energy injected by intense star formation. Here we discuss the possibility that the high 
surface brightnesses observed are not due to star formation, but due to cosmological gas accretion. We assume that all of the accretion 
energy is dissipated as shocks from the accreting gas. We show that in order to explain the high surface brightnesses both the mass 
accretion rate and energy would have to be much higher than expected from simulations or from equating the star formation with 
the accretion rate. We also investigate scaling relations between the surface brightness expected from accretion and for star formation 
through mechanical heating and photo-ionization, tiying to identify a regime where such accretion may become evident in galaxies. 
Unfortunately, the surface brightness necessary to detect the gas in optical line emission is about an order of magnitude lower than 
what has currently been achieved with near-infrared observations of distant galaxies. 
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1. Introduction 

Various processes have been proposed to explain how galaxies 
obtain the amount of gas necessary to fuel star formation and 
thei r growth. In the early hierarchical models of galaxy growth 
(e.g. lWhite & ReesI 19781) it was a combination of minor and ma- 
jor mergers, and gas cooling from a virialized dark matter halo, 
but it is now thought that much of the cooling may take place 
outside of the individual halos that subsequently accrete the gas, 
and that this gas will never be heated up to the virial temperature 
in the halo - such accretion mechanisms have been dubb ed "cold 
flows" (e.g. lDekel et alj|2009l;lva"n de Voort et al]|201 Ih . 

All these ideas rely on the basic assumption that baryons are 
strongly coupled to the dark matter through gravity and thus 
robustly follow the collapse of the dark matter. Without such 
coupling and the relative amounts of heating and cooling that 
this entails, neither cooling virialized halo gas nor gas from cold 
flows would deposit as much material over the same time range 
as they have been hypothesized to do. However, there is scant 
observati onal evidence for such cold flows or even hot halos. 
Recentlv. lAnderson & Brepmanl(l2010l) concluded that the miss- 
ing baryons are not in the halos of nearby galaxies and likely 
never were part of the Milky Way. In addition, there is evidence 
that gas in halos of distant gala xies is due to starb urst driven 
outflows and not to gas accretion ( ISteidel et al.ll20ldl) . 

Neverthe l ess, many mode ls (e.g. iDekel et all 120091: 

iBrooks et all l2009t iKeres et alJ l2009t) suggest that cooling 
baryons in the halos and/or cold accretion flows are needed 



* Data obtained as part of Programmes ID 074.A-9011, 075. A- 
0318, 075.A-0466, 076.A-0464, 076.A-0527, 076.B-0259, 077.B-0079, 
077.B-0511, 078.A-0055, 078.A-0600, 079.A-0341 and 079.B-0430 at 
the ESO-VLT. 



to explain the observations of distant galaxies, especially their 
large emission line wi dths, the apparent p revalence of rotating 
disks (as claimed by iForster Schreiber et a l. 2 006], and their 
high s tar-formation rates. On the other hand. Lehnert et aP 
(l2009h (hereafter L09) proposed that the large line widths 
observed are not due to gas accretion but to the intense star 
formation in the observed high redshift (z~2) objects. The 
same conclusion has subsequently b een reached by lOreen et al.l 
(l2010h . lElmegreen & Burkerd (l2010b (hereafter EB 10) modeled 
how accretion flows might explain the large line widths by 
dissipating the accretion energy as turbulence. They conclude 
that while large line widths are possible due to accretion, the 
phase when line widths are dominated by accretion must be 
short lived, about one dynamical time (~100 Myr), suggesting 
that intense, or more efficient, star formation may be more likely 
responsible for the large line widths. 



Despite this work, many subsequent studies have continued 
to suggest that gas accretion is a viable mechanism for explain- 
ing some of the properties of high redshift galaxies. In this pa- 
per, we examine the two following simple questions: (1.) Can 
the observed high Ha surface brightnesses in distant galaxies be 
explained by cosmological gas accretion at a rate sufficient to 
fuel their star formation? (2.) For what ratio of the star forma- 
tion and accretion rates would we expect the emission line lu- 
minosity generated from accretion to dominate? The impact of 
gas accretion can only be observed if the heating rate due to star 
formation is below some threshold. This threshold depends on 
the efficiency of (1.) the transformation of accretion energy into 
turbulence, (2.) gas heating and (3.) the transformation of the 
mechanical and ionizing energy of massive stars into turbulent 
and bulk motions, and heating. 
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Throughout the paper we adopt a flat Ho =70 km s ' Mpc ^ 
concordance cosmology with Qa = 0.7 and Q.m = 0.3. 

2. Ha luminosity and velocity dispersion 

We use ESO archival data from a variety of programs with 
SINFONI on the ESO-VLT of a sample of more than fifty galax- 
ies in the redshift range 1 .3-2.7, which will be described in more 
detail in Le Tiran et al. (2011, in preparation). The key to the 
sample is that all have sufficiently high Ha surface brightness 
to obtain spatially resolved line maps in a few hours integration 
time. The targets have Ho- surface brightness levels exceeding 
a few X 10"'** erg cm"^ s"' arcsec"' over more than 1-2 arcsec^ 
(10-20 kpc) in radius and total Ha luminosities above about 10^^ 
erg s"' . Due to cosmological dimming, there is a range of a fac- 
tor of ~5 in the faintest surface brightness line emission lev- 
els that can be probed as a function of redshift in our sample. 
However, it is worth repeating the conclusion from L09 that even 
our lowest redshift SINFONI sources are extreme compared to 
galaxies in the local volume, where only a few of the most in- 
tense starbursts have such high Ha surface brightnesses, and this 
on smaller physical scales only. 

Many studies which have analyzed the emission-line char- 
acteristics of distant galax ies examined velocity dispersions as 



an integrated property (e.g.lKhochfar & Sil k 2009; Burkert et al.l 
l2010t ICeverino et al. l2010l) . These analvses have tended to fa- 



vor gravity as the source of the observed characteristics, whether 
driven by gas accretion energy derived from the potential energy 
of the halo or from gravitational contraction of a gaseous and 
stellar disk. L09, on the other hand, focused on the relationship 
they found between spatially resolved Ha surface brightness 
and velocity dispersion in 10 high redshift galaxies observed 
with SINFONI, a subset of the sample of over 50 galaxies dis- 
cussed here. They suggest that neither smooth cosmological gas 
accretion, velocity dispersions driven by Jeans-unstable clumps, 
or turbulence generated by energy extracted by a collapsing disk 
can explain this relationship. They suggested instead that the 
velocity dispersion due to mechanical energy released by the 
observed intense star formation is of the form cr={e'EsFRy^^, 
where e is the efficiency at which the mechanical energy from 
a star-formation intensity is converted into turbulence and 
bulk flows in the interstellar medium of galaxies. This model has 
the advantage of having no free parameters since the efficiency 
can be constrain ed by observations of nearby galaxies or models 
et al.ll2006l and references therein). 
Rather than using spatiall y resolved spectrog raphic data for 
a pixel-by-pixel comparison. iGreen et al.l (1201 Ol) constructed a 
single, integrated Ha spectrum per galaxy, through a flux- 
weighted averaging of their pixels. Figure [1] shows the re- 
lationship between the total Ha luminosity and the mean 
velocity dispersion measured using integrated spectra derived 
with flux weighting per pixel, for our sa mple of high r edshif t 
galaxies, as well as for the galaxies from iGreen et al. (2010]), 
i.e., t heir selection of loc al SDSS galaxies. M51 (ITu11vu1974) . 
M82 (iLehnert et al.lll99 9 |) and diverse s a mples of h i gh red shift 
galaxies taken from'Epinat et al.J (l2009l) . iLaw et al.l (l2009h and 
[Lemoine-Busserolle et al. (2010). The relationship shows a sim- 
ilar trend as the one found between surface brightness and dis- 
persion per pixel in L09. For our purpose here, this remarkable 
similarity is important as it supports the notion that any hypoth- 
esis that explains the trend of increasing local optical emission 
line widths with local surface brightness, such as star formation 
(Lehnert et al. 2009), apparently also explains the similar rela- 
tionship between the integrated quantities (something that was 
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Fig. 1. Intrinsic mean Ha velocity dispersion (in kms ') as 
function of total Ha luminosity (in erg s"') for different sam- 
ples of local and high redshift galaxies (see § |2]i. The green 
line represents the Ha p roduced by a gas accretion model 
(iDekel. & Bimboimll2008[ see §|3]l. 



not pointed out by either L09 or lGreen et al 1120101) . We will dis- 
cuss this in greater depth in Le Tiran et al. (201 1, in prep.). 

3. Can cosmological gas accretion shocks produce 
the observed Ha surface brightnesses? 

L09 evaluated various mechanisms that can contribute to the 
large observed Ha line widths, such as clumps, gravitational 
instabilities, effects from the turbulence and smooth accretion, 
and suggested that only the mechanical output from star forma- 
tion could power such large random motions. This paper further 
examines the case for cold gas accretion and its observability in 
the z~2 Universe (e.g.;Keres et al. 2009; Dekel et al. 2009). Gas 
accretion onto galaxies is undoubtedly a complex process: this 
accreting gas would further cool through radiation and by gen- 
erating shocks upon hitting the ISM of the galaxy. The highly 
radiative shocks are an effective mechanism through which in- 
falling gas can lose its energy and momentum. But can we see 
direct evidence for this in our integral field Ha data sets? 

The underlying assumption that the infalling gas dissipates 
energy by colliding with gas already in a galaxy disk and that all 
the accretion energy, extracted from the halo potential, is con- 
verted into supersonic turbulence which then dissipates as fast 
shocks is undoubtedly an over-simplification, as in reality only 
part of the energy will go into turbulence and some of it into 
bulk motions, and shock velocities will span a wide range of 
velocities which is dependent on the structure (density and tem- 
perature) of the ISM. Although some, and perhaps most, of the 
energy will be dissipated in slower shocks in denser gas, the Ha 
luminosity of the shock will not be very high in shocks below a 
few tens of kms"' (ShuU & McKee 1979)." 

We can estimate the dissipation rate of the infalling gas us- 
ing a range of shock models, as we assume that the shock veloc- 
ities are equal to the observed velocity dispersions, which span 
the range of -50-250 kms ' in our sample. The sh ock mod- 
els are complementary in the sense that some (lAllen et al. 200^ 
focus on high velocity shocks where the precursor ionizing the 
gas in front of the shock is important (Vshock > 80-100 kms"'. 
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iDopita & Sutherland! 1995[). and others on lower velocity shocks 
( Ravmondlll979l:IShull & McKedll979h . 

This allows us to estimate the rate at which mass must pass 
through the shock in order to provide sufficient Ha surface 
brightness to be observable, as well as the total Ha luminosity 
produced by the accretion energy. Using the characteristics of 
the z~2 galaxies in our sample, we find a mean velocity disper- 
sion of almost 100 kms"' and a total Ho' luminosity of ~10"^^ '' 
erg s Correcting for the extinction typically observed for this 
kind of sources would increase their Ha luminosities by a fac- 
tor of 2-5 (L09). Using the shock models of Allen et al. (2008), 
we find that about 2% of the shock energy is emitted in Ha for 
a shock velocity of 100 kms independent of the gas density. 
The conversion rate of 2% is near the maximum efficiency for 
this process. 

An accretion energy extracted from the potential energy 
of the dark matter halo, of ~10'*-''^ erg s ' has been esti- 
mated for a baryonic accretion rate of about 100 Mq yr"' 
(iDekel. & Birn boim 2008). Our typical Ha luminosity, if en- 
tirely due to star formation, would be consistent with a star- 
formation rate of ~ 100 Mq yr"' for a Sal peter IMF with lo wer 
and upper mass limits of 0.1 and 100 Mq (lKennicutllll998L and 
a factor of a few lower for a more reasonable Kroupa IMF). The 
predicted Ha luminosity, if due to accretion, would then be about 
lO'*' '* (Macc/100 M0 yr"') erg s"', if the whole of the accretion 
energy is radiated through shocks with velocities of 100 kms 
with a conversion rate of 2%. Wit h 10% of the accretion energy 
going into supersonic turbulence (iKlessen & Hennebell3l2010h . 
gas accretion rates of many hundreds to thousands of M0 yr ' 
would be required for the total accretion energy to explain the 
Ha luminosity. Such high values for the gas accretion rate are 
well above what is commonly estimated from various models 
(e.g. iKeres et alj|2009t iGenel et aUlIOOSh . We demonstrate this 
in Figure [T] where we show the relationship between cr and Lfja 
for a gas accretion model. For this model, we have used the gas 
accretion rate as a function of virial velocity and the energy de- 
position rate due to accretion (Dekel, & Birnboim 2008), a con- 
version rate of 2%, and assumed the virial velocity is equal to 
the shock velocity. 

If gas accretion is truly happening at rates of up to several 
hundreds of M0 yr"' at z~2, it is likely that the integrated mea- 
sures of Ha luminosity and velocity dispersion are due to a com- 
bination of star formation and gas accretion. We can use the spa- 
tially resolved data to investigate if part of the surface brightness 
of these galaxies may be explicable by gas accretion. For this we 
will estimate the amount of mass that would need to be shocked 
to explain the observed surface brightnesses as a function of dis- 
persion. In Figure |2] we show the results for 10 representative 
galaxies in two redshift bins which, because of the strong de- 
pendence of surface brightness with redshift, allows us to show 
the full range of rest-frame Ha surface brightnesses within our 
sample of about 50 galaxies. We find that the surface bright- 
nesses can be explained by shock models, but only if the mass 
flow rates through the shocks are very high, with accretion rate 
densities of over a few hundred M0 yr ' kpc - implied for the 
highest surface brightnesses, whereas only for the lowest surface 
brightnesses a more reasonable level of 25-50 M0 yr"' kpc"^ is 
indicated. In all cases, to explain the surface brightness requires 
gas of reasonably high density (>100 cm"^). 

Taking a rate at which gas is shocked of 300 M0 yr"' kpc"^ 
would suggest 10'' M0 yr"' of sh ock heated gas, or ~ lO"' times 
the reasonable gas accretion rate (iDekel et al.ll2009b . To support 
star formation and to form a disk galaxy, the energy from ac- 
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Fig. 2. Observed pixel-by-pixel Ha line widths (assumed equal 
to the shock velocities) as function of Ha surface brightnesses 
(corrected for cosmological surface brightness dimming) in a 
sub-sample of 5 representative galaxies in each of two redshift 
bins, z~1.4 (+ symbols) and z~2.3 (x). Lines corresp ond to Ha 
surface brightnesses derived using shock models from lRavmondl 
( 1979 ) (black) and Allen et al. (2008) (blue). The numbers along 
each line indicate the surface accretion rates (in M0 yr"' kpc"^) 
which produce the Ha surface brightness. The 3 sets of lines are 
for pre-shock densities of 10, 100 and 1000 cm"^, respectively, 
from left to right. 

cretion must not be dissipated too rapidly least the gas loses too 
much angular momentum. Over a dynamical time of about 200- 
300 Myr (EE 10), the amount of gas that needs to be shock heated 
is of the order of 10'^^ M0, two orders of magnitude larger than 
the total ISM mass in a disk galaxy. 

4. Discussion and Conclusions 

The reason why it is likely difficult to note the effects of gas 
accretion on the ISM of presently observed galaxies is the rela- 
tively low efficiency of the conversion of accretion energy into 
line emission. An interesting question is under what circum- 
stances a galaxy might be observed that is undergoing a (for- 
mative) phase of intense gas accretion. 

We have attempted to quantify how the mass accretion rate 
would compare to the star-formation rate in both controlling the 
ionization of the gas and in mechanically exciting the gas (see 
Fig.O. We estimated that the Ha luminosity is pro duced at the 
rate of 2.6x10^" erg s"' Mq' ( iLeitherer etal.ll 19991 appropriate 
for a Kroupa IMF), that there is a 90% contribution of ioniz- 
ing photons (fy.sf), which implies an escape fraction of 10%, 
that accretion shock energy is converted into Ha luminosity at 
the rate of 2% {ef^""'''), and that 90% (face) of the accretion en- 
ergy is lost due to shocks which emit Ha at an efficiency of 
fflcc"^*^*- To compare this with the mechanical energy of the star 
formation we used a conversion factor of 7.4x10^*' erg s"' Mq 
dLeitherer et al.l[T999l) and a conversion efficiency of either 10 or 
20% (esp). 

This allows us to compare the relative efficiency of heating 
and ionization from the formation of massive stars with that of 
the mechanical heating through gas accretion. Can we identify 
areas in a plot comparing the rates of star formation and cosmo- 
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Fig. 3. Comparison of the star formation rates and cosmological 
gas mass accretion rates (in M0 yr"') at which accretion energy 
generates the same Ha luminosity as young stars through their 
ionizing photons (red solid and dashed lines) and as mechanical 
energy released by supemovae and stellar winds (blue lines); the 
difference between the two sets of lines is the efficiency at which 
the ionization or mechanical energy output couples to the ISM 
(see text for details). Both comparisons assume the same accre- 
tion model. Regions above the blue and red lines are dominated 
by the influence of the intense star formation. A comparison with 
the models by EBIO is also shown for 0.5, 1, 2, 3, 4 and 5 x 90 
Myr (from bottom to top, as indicated by the black arrow). The 
solid black line shows the situation when star formation and cos- 
mological gas mass accretion rates are equal (which the model 
reaches after a couple of dynamical times). 



logical gas accretion where observations may be made to dis- 
cover the duration and impact of gas accretion? To this end, we 
have used the simple model of the relationship between gas ac- 
cretion and star formation of EBIO. This model, which assumes 
a phase with constant accretion rate, allows us to compare, as 
a function of time, the expected total recombination line lumi- 
nosity due to infalling gas dissipating its energy through shocks 
versus that produced by the ionizing radiation and mechanical 
energy from star formation, whose rate is estimated by the total 
gas mass multiplied by the gravitational instability growth rate 
(see EBIO for details). 

In this comparison with the models of EBIO we see (Fig.[3]l 
that very quickly, within 45-90 Myr depending on the gas mass 
accretion rate, photo-ionization from massive stars dominates 
over the heating due to gas accretion. This implies that shock- 
like line ratios are not expected to be observed even if the en- 
ergy of the infalling gas dominates the dynamics of the ISM. 
After about 180 Myr, the mechanical energy of massive stars 
through supernovae and stellar winds begins to dominate, some- 
thing already concluded by L09. While we should not take the 
numbers literally, the analysis suggests that the time over which 
the effects of gas accretion will be evident is relatively short - 
less than a dynamical (orbital) time within a galaxy. 

So we expect to see the impact of gas accretion on either 
the ionization or dynamical state of the gas only early in the 
gas accretion process and only when the star-formation rates are 
relatively low. Within the context of the EBIO model, the star 
formation will roughly reach equilibrium with the accretion rate 



after a couple of dynamical times. If this were the case for our 
intensely star-forming galaxies at high redshift, we should not 
see any impact on either the surface brightness or dynamics of 
the emission line gas, as these are both controlled by star forma- 
tion, in agreement with wha t was first proposed by L09 and later 
supported by the findings of iGreen et al.l (l20I0h . 

Can we actually observe this early phase of gas accretion and 
its impact on either the ionization or the mechanical energy? The 
answer is likely no, at least not through Ha observations. Star- 
forming gala xies at high redshift have star- formation rates above 
10 Mq yr ' dForster Schreiber et"aDl2006l) . Comparing this with 
Figure [3] suggests there is only a small probability that the ef- 
fect of gas accretion can be observed. Taking a relatively op- 
timistic star-formation rate of less than 30 Mq yr"', which is 
equivalent to the Ha surface brightness generated by a gas ac- 
cretion rate of about 100 Mq yr using our scaling for mechan- 
ical energy, suggests that the average surface brightness over 1 
arcsec^ (barely enough to resolve a galaxy) would be about 4 x 
10"'^ erg s"' arcsec"- at z-2.2. This is about an order of mag- 
nitude less than what can be detected in a few hours observation 
with a spectrograph like SINFONI on an 8-lOm class telescope 
like the ESO-VLT, an d even beyond the longest integrations 
made with SINFONI dLehnert et all 1201 Ol) . Perhaps more un- 
fortunately, given the difficulty in designing spectrographs with 
large pixels on the next generation of large telescopes, reaching 
surface brightness levels beyond the current limits will remain 
difficult in the foreseeable future. 

Therefore, we conclude that currently, the impact of gas ac- 
cretion on galaxies at high redshift cannot be observed using the 
optical emission line gas. This is true especially for the early 
phases of galaxy formation where the gas accretion rate, despite 
it being thought of or der 100-300 M(7) yr ' for massive ga laxies 
and at high redshifts (iDekel et al.ll2009l: iKeres et alj|2009l) . will 
not be observable both because it does not liberate enough me- 
chanical energy to be observable with integral field spectrom- 
eters and because it is relatively quickly over- whelmed by the 
heating due to star formation, even at relatively modest rates. 
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